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ABSTRACT 


A  model  has  been  developed  for  the  prediction  of  the  time  dependent 
properties  of  the  gas  surrounding  a  water  droplet  under  irradiation  by  a 
high  intensity  laser  beam.  The  basic  phenomena  included  in  the  model  are 
absorption,  heating  and  vaporization  in  the  droplet  and  convection,  thermal 
conduction  and  diffusion  in  the  gas.  The  model  assumes  volume  absorption 
by  the  droplet,  a  constant  pressure  gas,  and  that  the  initial  vapor  velocity 
is  much  less  than  the  local  sound  speed.  The  conditions  under  which  these 
assumptions  are  appropriate  have  been  identified.  A  number  of  preliminary 
calculations  are  presented. 
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AEROSOL  PROPAGATION  EFFECTS 


I.  INTRODUCTION 

The  Aerosol  Propagation  Effects  program  is  concerned  with  the 
prediction  of  laser  beam  degradation  resulting  from  the  interaction  between 
a  high  power  laser  beam  and  atmospheric  aerosols.  The  particular  emphasis 
of  this  program  is  on  the  absorption,  scattering  and  thermal  blooming  of  a 
DF  laser  propagating  through  an  atmosphere  contaminated  by  hygroscopic 
aerosols.  The  overall  program  includes  a  review  of  the  properties  of 
atmospheric  aerosols,  a  study  of  the  atmospheric  heating  resulting  from 
acre  sol  absorption  of  laser  energy  and  first  order  predictions  of  the  laser 
bean  thermal  blooming  resulting  from  this  atmospheric  heating. 

The  particular  problem  addressed  in  this  mid-program  progress 
report  is  the  prediction  of  the  space/time  variation  in  the  temperature  and 
water  vapor  distributions  about  a  pure  water  droplet  irradiated  by  a  laser  beam. 
This  problem  is  relevant  to  beam  propagation  through  hygroscopic  aerosols 
on  humid  days.  Effects  resulting  from  the  varying  properties  of  real  aerosol 
solutes  and  original  Aitken  nuclei  will  be  discussed  in  the  final  report. 

The  description  of  the  general  problem  and  derivation  of  the  governing 
partial  differential  equations  for  the  droplet  and  surrounding  gas  will  be  given 
in  Section  II.  Specific  "steady  state"  and  time  dependent  computer  solutions  of 
the  problem  will  be  presented  in  Section  III.  The  conclusions  and  future 
directions  of  the  study  are  presented  in  Section  IV. 


II.  MODELING 

In  the  problem  under  consideration  the  aerosol  droplet  absorbs 
energy  and  is  heated.  As  the  droplet  temperature  rises  above  ambient 
the  droplet  conducts  heat  energy  to  the  surrounding  air  and  undergoes  surface 
vaporization.  The  vapor  leaves  the  droplet  with  a  directed  velocity  and  a 
pressure  wave  is  induced  in  the  surrounding  gas  in  response  to  the  mass  and 
heat  addition.  (If  the  rate  of  energy  absorption  were  sufficiently  large  the 
droplet  would  "explore",  creating  a  shock  wave;  however,  the  present  analysis 
will  be  limited  to  laser  intensities  and  droplet  sizes  such  that  the  directed 
velocity  of  the  vapor  leaving  the  droplet  surface  is  much  less  than  the  local 
sound  speed.  )  An  additional  feature  which  must  be  included  in  the  analysis  is 
the  fact  that  the  droplet  radius  decreases  with  time  due  to  mass  loss  resulting 
from  vaporization. 

A.  Gas  Equations 

A  number  of  approximations  have  been  made  to  simplify  the  solution 
of  this  problem.  The  most  prominent  of  these  is  that  the  gas  is  taken  to  be  a 
constant  pressure.  As  pointed  out  above,  pressure  gradients  are  induced  in 
the  gas  because  of  the  mass  and  heat  addition  from  the  droplet.  However,  if 
the  pressure  equilibrates  on  a  time  scale  fast  relative  to  the  times  for  conduction 
and  diffusion  then  the  heat  and  mass  transfer  may  be  assumed  to  be  occurring 
in  a  constant  pressure  gas.  The  characteristic  time  for  water  vapor  to  diffuse 
a  distance  L  in  air  is  approximately 

td  *»  l2/d  (1) 

where  D  is  the  diffusion  coefficient.  The  characteristic  time  for  heat  conduction 
in  air  is 

Tc  “  L'p  Cp/kA  <2) 

A 
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jvhere  P  is  the  density,  Cp  the  specific  heat  and  kA  the  thermal  conductivity 
of  air.  Lastly,  the  characteristic  time  for  pressure  equilibrium  is 


where  c  is  the  sound  speed  in  air.  Evaluating  these  expressions  with 
D  =  0.  24  cm2/sec,  PA  =  1.  2  x  10  3  gms/cc,  Cp,  =  0.  24  cal.  /gm-  K, 
kA  =  6.  2  x  10'5  cal.  / cm- sec -°K  and  c  *  3  x  1(T  cm/sec  results  in 


T 

D 


»  T 


P 


(4) 


for  L>  1  Mm.  The  characteristic  droplet  sizes  of  interest  in  this  work 
are  a  few  tenths  Mm  to  50  Urn,  and  for  the  typical  atmospheric  aerosol 
distributions  of  interest  the  average  spacing  between  droplets,  which  would 
be  the  heating  range  of  interest,  is  approximately  one  mm.  Thus,  it  would 
appear  that  the  approximation  of  a  constant  pressure  gas  is  reasonable. 

Another  simplifying  approximation  used  is  that  the  gas  properties  are 
taken  to  be  those  of  clean  air.  Although  this  approximation  would  suggest  that 
the  defining  equations  would  not  be  valid  for  high  water  vapor  concentrations 
this  is  not  the  case,  as  will  be  demonstrated  later  in  this  section.  Also,  the 
coefficient  of  thermal  conductivity,  kA  and  the  product  of  the  coefficient  of 
diffusion  and  number  density,  DN,  were  taken  to  be  constant.  These  quantities 
vary  by  <  15%  over  the  temperature  range  of  interest,  0  -  100°  C. 

The  last  simplifications  concern  the  fluid  mechan’.cs.  The  processes 
considered  in  the  analysis  are  thermal  conduction,  diffusion  and  convection. 
The  phenomena  of  thermal  diffusion  (Soret  effect)  and  thermo-diffusion 
(Dufour  effect)  have  not  been  included.  Generally,  these  processes  are  insig¬ 
nificant  in  flow  fields  of  the  type  under  consideration;  nonetheless,  the  validity 
of  this  assertion  will  be  determined  from  the  computer  results.  Furthermore, 
consideration  must  be  given  to  gravitational  effects.  As  the  gas  surrounding 
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the  droplet  is  heated  it  will  rise.  The  characteristic  time  for  the  gas  to  rise 
a  distance  L  is: 


Vg  -  {lUg  [ap/p])1/2 

2 

where  g  is  the  acceleration  due  to  gravity,  980  cm/sec  ,  and  Ap/p  is  the 

fractional  difference  in  density  between  the  heated  gas  and  the  surrounding 

medium.  A  strong  lower  bound  or  T  corresponds  to  Ap/p  ~  0.3.  Even 

8 

in  this  limit  it  can  be  shown  that 


(5) 


T 
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T  « 


c 


T 
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for  L  less  than  one  mm.  Of  course  this  effect  might  have  to  be  included  for 
long  time  (t  >  0.  1  sec)  irradiation.  However,  in  this  case  it  would  perhaps 
be  more  appropriate  to  treat  the  aerosols  as  an  aggregate  rather  than  individually. 

Within  the  framework  of  the  above-mentioned  approximations,  the  general 
relationships  describing  the  system  are 


a)  Conservation  of  mas  s 


A  p  /A  t 


-  V  •  (P  V) 


(6) 


where  V  is  the  gas  velocity 

b)  Euler's  equation 


A 


V/A  t 


(7) 


c)  Conservation  of  energy 


A  (1/2  p  V  +  pC  T)/A  t  = 


P  V  (1/2  V  + 


CpT)  - 


7  T 


(8) 


where  Cv  is  the  specific  heat  at  constant  volume 


and 


d)  The  diffusion  equation 

6  P  /£>t  =  v  •  (D?  p  -  p  V)  (9) 

where  p^  is  the  water  vapor  density  and  D  is  the  H^O-air  diffusion 
coefficient. 


If  the  droplet  is  heated  uniformly  (see  Section  B),  the  system  is 
spherically  symmetric  and  it  can  be  shown  that  the  ordered  kinetic  energy 
terms  drop  out  of  the  equations  under  the  assumption  of  constant  pressure. 
The  relevant  gas  equations  for  r  >  R,  where  R  is  the  droplet  radius,  are 


bT/bt 


-  V  bT/br  + 


k  k  A  T 
A _ 

r2Cp  PM 
A 


A 


£  (r2Wbr), 


(10) 


TT  (r2  Xv/6r),  (11) 

b  2 

—  (r  bT/br  ),  (12) 

where  k  is  Boltzmann's  constant,  A  is  Avogadro's  number,  P  is  the  total 
pressure  taken  as  one  atm.  ,  M  is  the  molecular  weight  of  air,  N  is  the 

A 

total  number  densi.y,  and  X  is  the  mole  fraction  of  water  vapor. 

These  equations  have  been  formulated  in  terms  of  number  density  rather 
than  mass  density,  since  in  a  constant  pressure  flow  of  two  gases  with  different 
molecular  weights  the  quantity  NT,  rather  than  P  T,  remains  constant.  As 


bX  /bt  =  -  VbX  /br  +  -- 

V  V  Pr2 


b  V/b  r  =  -  2  V/r  + 


k  k  A 
A 


r 
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mentioned  earlier,  these  equations  were  derived  under  the  assumption 

X  «  1.  In  reality,  the  quantity  Cp  M  should  be  replaced  by  the 
v  A  A 

expression 

C  MA  (1  -  X  )  +  C_  M  X 
P  A  v  P  v  v 

A  v 

where  the  subscript  v  refers  to  water  vapor,  and  the  diffusion  coefficient  should 

be  replaced  by  that  for  full  binary  diffusion.  However,  over  the  temperature 

and  water  vapor  ranges  of  interest  in  this  problem  these  latter  two  quantities 

are  reasonably  constant  and  well  represented  by  Cp  and  D,  respectively. 

A 

Thus  Eqs.  (10)- (12)  are  approximately  valid  for  water  vapor  concentrations 
approaching  unity. 

B.  Droplet  Equations 

In  the  present  analysis  the  temperature  field  across  the  droplet  has 
been  taken  to  be  uniform.  As  long  as  the  laser  radiation  is  absorbed  volumetrically 
by  the  droplet,  this  approximation  is  reasonable  until  the  droplet  temperature 
approaches  the  boiling  point  of  water  (see  Section  III).  The  absorption  co¬ 
efficient  of  a  water  droplet  is  a  complex  function  of  droplet  radius,  laser 
wavelength  and  the  complex  index  of  refraction  (see  for  example  Ref.  1). 

2 

However,  for  4  lim  radiation  liquid  water  has  a  bulk  absorption  coefficient, 
a,  of  -  100  cm'1,  and  it  can  be  shown  that  for  these  conditions  water  droplets 
of  radius  <  50  Um  undergo  volume  absorption.  (The  general  relationship  for 
the  absorption  coefficient  of  a  water  droplet  may  be  found  in  Ref.  3  and 
references  therein.  ) 

The  basic  equatic  n  connecting  the  droplet  and  gas  dynamics  is  conserva¬ 
tion  of  energy  in  the  droplet. 
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4/3HR  la  =  4/3  TT  R  P 


C  d  T_ /d  t 
D  PD  D 


+  4TTR2kA(-6T/0r|R)  ^ 

,  M  I 

+  4  IT  R  AH  — ^  (N  V)  _ 

v  A 


The  LHS  of  Eq.  (13)  is  the  rate  of  absorption  of  energy,  where  I  is  the  laser 
intensity,  and  the  RHf  represents  the  energy  disposition  into  droplet  heating, 
conduction  and  vaporization  respectively.  The  quantity  AH^  is  the  change  in 
enthalpy  required  to  proceed  from  liquid  to  gaseous  state  and  the  subscript  D 
refers  to  droplet  properties.  Note  that  the  temperature  is  continuous  across 
the  droplet  surface;  i.  e.  TD  =  T  (R).  The  terms  involving  the  ordered  kinetic 
energy  of  the  vapor  are  negligible  for  the  present  considerations  and  have  been 
left  out  of  Eq.  (13). 

The  vapor  concentration  at  the  surface  of  the  droplet  has  been  related  to 
the  droplet  temperature  by  the  Clausius-Clapey ron  relationship. 


d  In  X 
_ v 

d  T 


P 


A  H 


R  T 
v 


(14) 


where  R  is  the  ideal  gas  constant  per  unit  mass  of  water  vapor.  It  has  been 
v  . 

demonstrated4  that  this  relationship  properly  describes  the  surface  vapor  pressure 
of  a  water  droplet  when  the  vapor  velocity  is  less  than  the  sound  speed.  In  the 
present  analysis  Eq.  (14)  has  been  used  in  its  integral  form 

X  (R)  =  X  exp 
v  v_ 


(15) 


uhere  the  subscript  ®  refers  to  ambient  conditions.  Eq.  (IS)  is  derived  under 
the  assumption  that  AH  is  a  constant.  Indeed  A  H  varies  by  less  than  10^ 

V  V 

between  0-100  C  and  has  been  taken  to  have  the  constant  value  of  585cal/gm. 

The  last  equation  connecting  the  droplet  and  pas  behavior  is  the  con¬ 
servation  of  n  ass  of  the  vapor.  If  a  spherical  shell  of  thickness  A  R  is 
constructed  about  the  droplet,  then  it  is  required  that  the  rate  of  mass  loss  of 
the  droplet  be  balanced  by  the  rate  of  increase  of  vapor  within  the  shell  plus  the 
rate  of  vapor  leaving  the  exterior  of  the  shell.  With  some  straightforward 
mathematical  manipulation  this  leads  to  the  relationship 


(NV) 


R 


-  (1  -  X  (R)  ) 
v 


DN  A  X  It  r 
v 


(16) 


C .  Effect  of  Varying  Droplet  Size 

So  far  in  the  analysis  the  effect  of  decreasing  droplet  size  has  not  been 
considered.  This  effect  can  be  uncoupled  because  the  rate  of  change  of  droplet 
size  is  small  compared  to  the  rate  of  change  of  the  gas  parameters  (because  of  the 
large  difference  between  liquid  and  gas  densities).  Of  course,  it  has  been 
implicitly  assumed  that  the  time  varying,  rather  than  initial,  value  of  R  will 
be  used  in  the  solution  of  the  equations. 

The  governing  equations  have  been  recast  in  terms  of  the  variable 

Z  =  R  (t) / r  (17) 


In  terms  of  this  variable  the  equations  may  be  written  as 


6  T/ist  =  R‘  1  (QZ4  -  Z  dR/dt)  iT/iZ 


kkAA  -242  ,„2 

cT  PM  TR  z  &  Tltz 
PA  A 


(18) 
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6  X  /6t  =  R"  1  (QZ4  -  Z  dR/dt)  b  X  /b  Z 

v  v 


♦  Tr'V^X  /^.2 

P  V 


(19) 


fi  Q/fiZ 


k  kA  A 

A 


- 1  .  2  JX  2 

R  6  T/6Z 


(20) 


wh  ere 


P  M 


dR/dt  =  - 


AODkT 


Q  (Z  =  1) 


(21) 


and 


Q  =  V/Z 

The  transformed  boundary  conditions  at  r  =  R  (Z  =  1)  are 


(22) 


la  =  o  _  C  6  T/bt 
D  '  D 


+  3  R"  k  6  T/&  7 
A 


Z=1 


Z=  1 


.MAH  P 

+  3R-  Q(Z-l). 


(23) 


Q  (Z  =  1 )  =  R 


1  k  T  DN 


r  1  - 1  i 

II 

N! 

X 

i 

&X  /6Z 

L  v  J 

V 

(24) 


Z=1 


and  the  Clausius-Clapeyron  relationship,  Eq.  (15). 


HI.  PRELIMINARY  CALCULATIONS 

No  analytic  solution  has  been  found  for  the  system  of  Eqs.  (IS), 

(18)  -  (24).  A  computer  program  has  been  developed  for  their  numerical 
evaluation  which  employs  a  fully  implicit  Crank  -  Nicolson  finite  diffen  r  ci 
scheme.  A  complication  in  using  me  /.  co-ordinate  is  that  a  constant  grid 
size  in  /.  space  results  in  a  limited  number  of  grid  points  at  large  r,  where 
the  solution  may  be  most  interesting.  This  complication  has  been  circum¬ 
vented  by  modifying  the  finite  difference  scheme  to  allow  for  a  change  in 
grid  spacing  at  small  /. .  In  the  present  mode  of  operation  the  grid  spacing 
tor  0.  0  <  /  <  0.  05  is  an  order  of  magnitude  finer  than  that  for  0.  05  <  / 

L  1.00. 

The  computer  code  is  now  operational  and  a  number  of  preliminary 
calculations  have  been  performed.  However,  those  are  insufficient  in  a 
number  to  provide  an  overview  of  the  problem.  A  "steady  state"  analysis 
of  the  equations  lias  simultaneously  been  performed  in  order  to  predict 
the  effects  of  varying  fundamental  parameters  such  as  particle  radius  and 

laser  intensity. 

A.  Steady  State  Solutions 

In  the  limit  where  the  velocity  term  in  Eqs.  (18)  and  (19)  can  be 
neglected,  the  equations  effectively  uncouple,  and  it  can  be  readily  demon¬ 
strated  that  the  c»eady  state  solutions  for  A  T(  T  -  T®  )  and 
^  X  (  X  -  X  )  scale  inversely  with  r,  i.  c.  , 

V  V  V  00 

A  T( r )  =  (R/r)  AT  (R)  (25) 

and 

A  X  (r)  (R/r)  AX  (R)  (2b) 

v  v 

Indeed,  even  the  time  dependent  solutions,  which  may  be  found  in  standard 
pie,  Ref.  5)  vary  approximately  as  r  for  r  ~  R.  Thus, 


texts  (for  exam 


in  the  limit  of  lew  velocity  the  gradients  app  aring  in  the  equation  for  drop¬ 
let  energy  conservation,  Eq.  (23),  may  be  evaluated  by  use  of  Eq.  (2C)  and 
(2b). 

If  this  is  done,  use  of  Eqs.  (15),  (23)  and  (24)  results  in 


d  (  A  TD)/dt  =  — j 


R  1  a  -  k.  AT 


R  C  D  CP 


D 


D 


AH  M 

v  v 

A 


DN(  1  -  X  (R)  ) 
v 


(X  (R) -  X 

v 


(27) 


where 


X  (R)  -  X 

V  v 


X 

Vcr, 


*Hv  *TP 

Rv<AV 


} 


Although  this  differential  equation  is  separable,  it  affords  no  simple  analytic 
solution.  One  feature  of  interest  is  the  maximum  value  of  ATp  which  is 
reached  when  *he  time  derivative  is  zero.  Equation  (27)  in  this  limit  re¬ 
duces  to  a  transcendental  equation  for  ATD_max  which  may  bc  evaluatc(1 
as  a  function  of  the  parameter  R^Ia  .  The  solution  to  this  equation  is  shown 
in  Fig.  1  over  a  range  of  5  orders  of  magnitude  in  the  parameter  R  I  a. 
Shown  for  comparison  are  the  predicted  maximum  droplet  temperatures 
for  the  related  cases  of  "heat  conduction  only"  and  "vaporization  only". 

These  calculations  were  performed  for  initial  conditions  of  T^,  293  K, 

-2 

X  =  2.  1  x  10 

v®  2 

Note  that  for  the  full  range  of  R  Ijl  shown  the  ATD_max  predicted 

for  the  case  of  heat  conduction  alone  is  sigrificantly  higher  than  that  pre¬ 
dicted  when  both  heat  conduction  and  vaporization  were  included.  On  the 
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other  hand,  the  predictions  for  the  case  of  "vaporization  cnl>  re  relatively 

close  to  the  full  solution.  The  reason  for  this  is  that  the  major  portion,  ol 

tht  absorbed  energy  goes  into  vaporization  rather  than  heat  conduction.  This 

is  shown  explicitly  in  Fig.  1  where  the  fraction  of  the  absorbed  energy  whicl 

2  2 

is  lost  to  heat  conduction  is  plotted  vs.  R  la.  For  low  R  Io  this  traction 

2 

asymptotes  to  ~  0.3,  and  falls  off  rapidly  with  increasing  R  la  .  This 
rapid  fall-off  occurs  as  the  droplet  temperature  approaches  the  boiling 
point.  This  is  because  the  vaporization  rate  increases  exponentially  in  this 
region,  whereas  the  rate  of  heat  conduction  approaches  a  constant.  The 
relative  value  of  this  fraction  is  determined  by  the  ratio  of  the  last  two 


terms  on  the  RHS  of  Eq.  (27). 


For  small  A  T_  this  ratio  reduces  to 

D-  max 


Power  into  heat  conduction 
Power  into  vapo  Tzation 


k  A  R  T  '  A  ( 1  -  X  ) 
A  v  ®  v  t> 


X  ,DN  (  A  H  )  M 

Vq3  V  V 


»  0.43 


(28) 


in  agreement  with  the  asymptotic  value  shown  in  Fig.  1  (i.  e.,  0.43/1. 43  »0.  •>) 
It  should  be  pointed  out  that  as  the  droplet  temperature  approaches  the 
boiling  point  the  convective  terms  left  out  of  the  steady  state  analysis  become 
more  important  and  convective  heat  transfer  might  dominate  conductive  heat 
transfer.  The  steady  state  equations  including  convection  are  more  compli¬ 
cated,  but  can  be  analyzed.  A  discussion  of  convective  effects  will  be  in¬ 
cluded  in  the  final  report. 

The  fact  that  vaporization  dominates  heat  conduction  as  an  energ 
loss  mechanism  allows  an  approximate  analytic  formulation  for  several  oi 
the  variables  of  interest  in  the  problem.  In  particular,  in  steady  state  it 
can  be  readily  shown  from  Eq.  (27)  that  the  rate  of  mass  loss  from  the 
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droplet  can  be  expressed  as 


dm/dt  =  -4/3  >rR3IaMH^  (29) 

or  in  terms  of  the  droplet  radius 

dR/dt  =  -RIo/(3AH  p  )  (30) 

From  Eq.  (30)  it  can  be  seen  that  the  characteristic  e-folding  time  for  drop¬ 
let  radius  is 

t  R  =  3AHv0p/Ia  (31) 

which  is  inversely  proportional  to  I  a  and  independent  of  initial  particle 
size.  Lastly,  the  vapor  velocity  at  the  droplet  surface  is  given  by 

V(R)  **  M  RIa/(3M  piH  )  (32) 

A  V  Vf 

Note  this  latter  quantity  scales  as  RI  a  rather  than  R  la,  and  thus  is  less 
sensitive  to  changes  in  radius  than  the  droplet  temperature. 

A  last  point  before  leaving  the  discussion  on  steady  state  solutions 
concerns  the  importance  of  temperature  gradients  within  the  droplet.  If 
one  considers  volume  heating  followed  by  heat  conduction  as  the  only  phen¬ 
omenon  occurring  within  the  droplet,  the  defining  steady  state  energy 
equation  becomes 

4/3TTr3Ia=-4TTr2kD&T/£>r  (r-R;  (33) 
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where  k  is  the  coefficient  of  thermal  conductivity  of  liquid  water,  taken 
D  j  ^ 

as  1.4  x  10"5  cal/cm-sec-°K.  Eq.  (33)  requires  that 


R2Ig 
6  k 


(30 


where  the  subscrips  C  and  S  refer  to  conditions  at  the  droplet  center  and 
surface,  respectively.  This  temperature  variation  is  also  shown  in  Fig.  1 
and  can  be  seen  to  be  significantly  lower  than  A  Tp  max  unt^  t*le  droplet 
approaches  the  boiling  point.  Of  course  Eq.  (34)  is  not  realistic  in  that 
region  since  it  would  require  the  droplet  to  superheat.  In  reality,  vapor¬ 
ization  would  occur  within  the  droplet  and  original  nuclei  effects,  etc.  could 
be  quite  important.  In  any  event,  it  is  clear  that  the  approximation  of  a 

constant  temperature  field  across  tin  droplet  is  not  approp’  iate  at  large 

„  2 

values  of  R  la. 

D.  Time  Dependent  Solutions 

The  full  finite  difference  scheme  has  been  exercised  for  a  limited 

number  of  cases.  As  an  example  of  the  results,  droplet  temperature  his- 

5 

tories  vs.  time  are  shown  in  Fig.  2  for  an  la  of  2.  5  x  10  ral/cc-sec 
and  two  initial  particle  sizes,  R  =  3,  10p.m.  Note  that  the  droilet  tempi  ra- 

tures  rapidly  rise  to  a  peak  value  and  then  slowly  decay.  This  decay  results 
from  the  decrease  in  droplet  radius  with  time.  The  radii  of  both  the  3  \nd 
10  pm  particles  decrease  by  ~  20 To  in  2  msec  and  exhibit  a  characteristic 
e-tolding  time  as  given  by  Eq.  (31).  The  calculated  droplet  temperatures 
from  their  peak  values  onward  are  in  good  agreement  with  the  steady  state 
predictions  shown  in  Fig.  1. 

Examples  of  the  gas  radial  temperature  distributions  for  several 

configurations  are  shown  in  Fig.  3.  Curves  C  and  D  are  for  the  same 

5  3 

case,  R  =  3am,  la  =  2.  5  x  10'  cal/cm  -sec,  at  two  different  times. 

Note  that  the  profiles  scale  approximately  inversely  with  r  until  large  r, 
where  they  drop  off  more  rapidly.  At  small  r  the  temperature  profile  is 
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t,  msec. 


2 


Fig.  2 


The  Rise  in  Droplet  Temperature  Over  Ambient  vs.  Time  for 
Two  Different  Droplet  Sizes. 
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1 

lower  at  the  longer  time.  This  is  because  of  the  droplet  cooling  which 
occurs  as  a  result  of  decreasing  droplet  size.  At  large  r  the  situation  is 
reversed  since  the  gas  temperature  in  that  region  is  still  being  driven  up¬ 
wards  towards  an  r’1  distribution.  The  characteristic  time  to  approach 
this  distribution  at  any  distance  r  is  given  approximately  by  Eq.  (2),  and 
is  independent  of  the  intensity  of  the  incident  irradiation. 

5 

Curve  B  is  a  similar  calculation  for  a  10  pm  particle,  la  =  2.5  x  10 
cal /cc- sec,  at  a  time  of  2  msec.  Once  again,  tho  profile  scales  approxi¬ 
mately  as  r’1  for  small  r.  Shown  for  comparison  as  curve  A  is  another 

6 

calculation  for  a  10  p  droplet  with  an  I  CL  of  2,  5  x  10  cal/cc-sco  at  a 

time  of  0.  2  msec.  Curves  A  and  B  correspond  to  two  cases  with  the  same 

energy  input  (IaAt  ).  However,  they  are  radically  different.  Curve  A  is 

sharply  bowed  away  from  an  r  1  distribution,  demonstrating  the  effect  of 

the  convection  term  in  Eq.  (18).  It  is  recalled  from  Eq.  (32)  that  the  initial 

vapor  velocity  scales  linearly  with  la  and  is  thus  an  o-der  of  magnitude 

larger  in  case  A  than  in  case  B.  On  the  other  hand,  curve  A  falls  below  B 

1/2 

at  large  r,  since  the  "conduction  length"  of  the  gas  scales  as  t  ,  Eq.  (2). 

No  calculations  of  the  radial  water  vapor  (AX  )  distributions  are 

v 

presented.  However,  these  are  quite  similar  in  shape  to  the  temperature 
distributions.  This  is  to  be  expected  since  it  was  demonstrated  earlier 
that  the  characteristic  time  for  diffusion  was  approximately  the  same  as 


that  for  heat  conduction. 


IV.  SUMMARY  AND  FUTURE  DIRECTIONS 

A  model  has  been  developed  for  the  prediction  of  the  time  dependent 
behavior  of  the  gas  surrounding  a  water  droplet  under  irradiation  by  a  high 
intensity  laser  beam.  The  basic  phenomena  included  in  the  model  are  ab¬ 
sorption,  heating  and  vaporization  in  the  droplet  and  convection,  thermal 
conduction  and  diffusion  in  the  gas.  The  model  assumes  volume  absorption 
by  the  droplet,  a  constant  pressure  gas,  and  that  the  initial  vapor  velocity 
is  much  less  than  the  local  sound  speed.  The  conditions  under  which  these 
assumptions  arc  appropriate  have  been  identified. 

A  computer  code  employing  a  iully  implicit  Crank- Nicolson  finite 
difference  scheme  has  been  developed  for  the  numerical  evaluation  of  the 
governing  differential  equations  and  a  number  of  preliminary  calculations 
has  been  performed.  Simultaneously,  a  steady  state  analysis  of  the  system 
has  been  performed  in  order  to  identify  the  salient  features  of  the  modeled 
phenomena  over  a  wide  range  of  the  problem  parameters.  From  this  ana¬ 
lysis  it  was  demonstrated  that  the  model  was  not  strictly  appropriate  at 
large  values  of  the  quantity  R2  I  a  ,  where  the  droplet  temperature  approaches 

the  boiling  point. 

A  number  of  tasks  are  presently  in  process  on  this  program  and  re¬ 
sults  will  be  included  in  the  final  report.  These  include: 

A)  A  review  of  the  properties  of  atmospheric  aerosols,  including 
particle  size  distribution  functions,  absorption  coefficients,  original 
nuclei,  etc. 

B)  An  extension  of  the  steady  state  analysis  to  include  convection 

effects . 

C)  A  matrix  of  computer  calculations  of  the  time  dependent  radial 
temperature  profiles  about  a  droplet.  Parameters  to  include  particle 
radius,  laser  intensity,  etc. 

D)  First  order  estimates  of  the  thermal  blooming  arising  in  typical 
scenarios. 
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E)  A  study  of  the  effects  of  original  Aitken  nuclei,  phase  changes 
associated  with  original  nuclei  material  and  production  of  new  chemical  species. 
Possible  internal  boiling  effects  with  accompanying  superheating,  fragmenta¬ 
tion,  mass  transport  and  shockheating  of  the  air  surrounding  the  droplet  vwll 
also  be  examined. 


-20- 


ACKNO  WLEDGEMENT 


The  authors  wish  to  acknowledge  useful  discussions  with 
Dr.  J.  D.  Teare  and  Mr.  P.  F.  Lewis.  Mr.  R.  Lundberg  was  responsible 
for  the  computer  programming. 

This  research  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by  ONR  under 
Contract  No.  N000  1 4- 74-C-0254. 


-21- 


REFERENCES 


1.  Van  de  Hulst,  H.  C.,  "Light  Scattering  by  Small  Particles ", 

John  Wiley  and  Sons,  New  York  1957. 

1.  Hale,  G.  M.  and  Qucrry,  M.  R.  ,  "Optical  Constants  of  Water 

in  the  200-nm  to  200-Lm  Wavelength  Region",  App.  Opt.  12, 

555  (1973) 

3.  Hodges,  J.  A.,  "Aerosol  Extinction  Contribution  to  Atmospheric 
Attenuation  in  Infrared  Wavelengths",  App.  Opt.  1  1 ,  2304  ( 1972). 

4.  Williams,  F.  A.  .  On  Vaporization  of  Mist  by  Radiation",  Int. 

J.  Heat  Mass  Transfer  8,  575  (1  965). 

5.  Carslaw,  H.  S.  and  Jaeger,  J.  C.  ,  "Conduction  of  Heat  in  Solids , 
Oxford  at  the  Clarendon  Press,  1959. 


-22- 


